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ABSTRACT 

We present Subaru spectroscopy of the early-type lensing galaxies in the lens systems 
HST 14113+5211 (at the redshift z L = 0.464) and B 2045+265 (z L = 0.868), being 
aimed at measuring the velocity dispersion of the lensing galaxies as an important 
component to their mass distributions and internal dynamics. For HST 14113+5211 
we have obtained 174 + 20 km s _1 (1 a) inside an aperture of 0."9, and for B 2045+265 
we have obtained 213 ± 23 km s _1 inside an aperture of l."2. To extract the signif- 
icance of these information on the mass distributions and stellar dynamics of the 
lensing galaxies, we construct two alternative models for both the mass and velocity 
distributions. It is found that the mass-to-light ratios as derived from a joint lens- 
ing and stellar dynamics analysis are virtually consistent with those obtained from 
the use of the Fundamental Plane. Also, for HST 14113+5211, we find that the total 
mass distribution is well reproduced by a power-law density profile with an index of 
2, thereby suggesting that a singular isothermal model is a good fit to the lensing 
galaxy. In contrast, the corresponding slope for B 2045+265 is flatter than isothermal, 
suggesting that additional contributions to lensing mass from surrounding structures 
(possibly the group of galaxies) provide the observed angular separations of multiple 
images. 
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1 INTRODUCTION 

Massive early type galaxies are predicted to form at the 
highest density peaks in the early universe and to evolve 
through subsequent mass assemblies. Measuring their dy- 
namical structure and stellar content over the cosmic time 
provides important constraints not only on their formation 
history itself but also on models of structure formation sce- 
nario and cosmic star formation history. 

Studying the internal structure of early type galaxies 
is, however, difficult because of lack of dynamical tracers 
at large radii such like Hi gas in spirals, and of degeneracy 
between kinematic properties of dynamical tracers and mass 
distributions: The stellar dynamics suffers from the mass- 
anisotropy problem, and the gravitational lensing does from 
the mass-density profile problem as well. 

Those degeneracies can be broken by combining those 
two probes, since they nicely complement each other (Treu & 
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Koopmans 2002; 2004 (hereafter TK04); Koopmans & Treu 
2003). A combined analysis thus places important limits on 
the distributions of luminous and dark matter in a lensing 
early type galaxy. 

We select two lens systems, HST14113+5211 (zl = 
0.464; Fischer, Schade & Barrientos 1998; Lubin et al. 2000) 
and B 2045+265 (z L = 0.868; Fassnacht et al. 1999), where 
the lensing galaxies are an early-type galaxy. We conducted 
the direct measurement of velocity dispersion of the lensing 
galaxies of those systems using the Faint Object Camera and 
Spectrograph (FOCAS, Kashikawa et al. 2002) mounted on 
the Subaru Telescope. We perform a joint gravitational lens- 
ing and stellar dynamics analysis to explore the luminous 
and dark matter distribution of lensing galaxies focusing 
on the slope of (dark and total matter) density profile and 
mass-to-light ratio. Also, we utilize the fundamental plane 
relation of early type galaxies for deriving an alternative es- 
timate of the mass-to-light ratio, which, after combined with 
the joint analysis, provides a useful constraint on the mass 
distribution of the lensing galaxy. 

The paper is organized as follows. In §2, we show the 
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observations and data reduction. In §3, we derive the mass- 
to-light ratio of the lensing galaxies from the fundamental 
plane relation of early type galaxies. In §4, models for grav- 
itational lensing and stellar dynamics are presented. §5 is 
devoted to the results of our model analysis, and discussion 
and concluding remarks are drawn in §6. 

Throughout this paper, we adopt fio = 0.3, Ao = 0.7, 
and h = H /W0 km s" 1 Mpc" 1 = 0.65 for the relevant 
estimations. 

2 VELOCITY DISPERSION MEASUREMENT 
2.1 Observations 

Spectroscopic observations were made with the Subaru 
8.2m telescope (lye et al. 2004). The FOCAS spectrograph 
(Kashikawa et al. 2002) was configured with a 0". 4- width 
slit, a 300 grooves mm -1 grism which gives 1.40 A per 
pixel, and a Y47 order-cut filter to obtain optical spec- 
tra of the lensing galaxies of gravitational lens systems 
HST 14113+5211 and B 2045+265. The on-chip binning was 
set to 3 (along spatial direction to give 0".3 per pixel) by 1 
(along wavelength direction). 

For HST 14113+5211, observation was made on 2002 
June 13 and 14 (UT). The slit was placed along the major 
axis of the lensing galaxy, whose position angles (PAs) was 
PA= 38.1°. We obtained ten 1800 seconds exposures, and 
the total exposure time was 5 hours. For B 2045+265, obser- 
vation was made on 2002 June 13 (UT). The slit was placed 
along the major axis of the lensing galaxy (PA= 120.1°). We 
obtained seven 1800 seconds exposures, and the total expo- 
sure time was 3.5 hours. The seeing conditions were around 
0.4" - 0.8" during the two observing nights. 



2.2 Data reduction 

Basic data reductions were made following Ohyama et al. 
(2002). In the followings, we focus on the spectra of the 
lensing galaxies. Wavelength accuracy and resolution of the 
galaxy spectra were measured, within the wavelength region 
for the Fourier cross-correlation analyses (see section 2.4), by 
means of Gaussian fittings of several narrowest (unblended) 
sky OH lines. Here all measured values are shown in the 
rest frame of each galaxy, since all the following velocity 
dispersion measurements were made in their rest frames. For 
the spectrum of HST 14113+5211, we found that wavelength 
resolution is 4.5A FWHM at 4000 - 4450A (the blue fitting 
region) and 3.2A at 4780 — 5170A (the red fitting region. 
See below for details of the spectrum fitting regions). The 
wavelength accuracy is found to be typically ~ 0.09A in 
RMS over 3800 - 6500A, and it gets slightly worse to ~ 
0.28A at the bluest wavelength (~ 4000A) for the cross- 
correlation analyses. For B 2045+265, wavelength resolution 
is 2.7A FWHM at 4100-4500A, and wavelength accuracy is 
typically ~ 0.07A in RMS over 3800 - 4500A. Aperture size 
to extract an 1-dimensional galaxy spectrum are 0".9 (3- 
binned pixel, corresponding to 8.3 kpc at z — 0.464 for the 
adopted cosmology) and l".2 (4-binned pixel, 19 kpc at z = 
0.868) for HST 14113+5211 and B 2045+265, respectively, 
to obtain spectra with maximum signal-to-noise quality. 
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Figure 1. Continuum-normalized spectrum of HST 14113+5211 
(black) in the rest frame of the galaxy. The continuum-normalized 
spectrum of a template star, HD 126778, after Gaussian- 
broadened to match the galaxy spectrum (<j = 174 km s — 1 ), 
is shown (blue) with a vertical offset of —0.3 for clarity of the fig- 
ure. Residual spectrum (galaxy — Gaussian-broadened template) 
(red) and a sky spectrum (black) are also shown. Here the sky 
spectrum is shown, with a vertical offset of —0.5, to identify wave- 
length at which strong sky emission could damage the observed 
galaxy spectrum, and is flux-scaled by an arbitrary number. At- 
mospheric extinction features are marked with vertical black bro- 
ken lines. Two wavelength regions used in the analyses of the ve- 
locity dispersion (blue and red fitting regions) are marked by two 
vertical green lines (see a main text for details of the wavelength 
fitting regions). Prominent features seen in the galaxy spectrum 
are identified as shown in upper portion of the figure. 

2.3 Basic properties of the spectra 

Reduced rest-frame spectra, after continuum normalization, 
are shown in Figures[T]and|21 For HST 14113+5211, deep Ca 
II H and K absorption lines and prominent G band feature 
are evidently seen in the spectrum, suggesting that major 
contribution in the observed wavelength region could be at- 
tributed to the late type giant stars (late G giant - early K 
giant stars). The redshift of the lensing galaxy is measured 
to be z = 0.4644+0.0002 from these prominent features, and 
is consistent with previous measurement (Lubin et al. 2000). 
In the redder part of the spectrum, several more features, in- 
cluding H/3 and some Fe absorptions (Fe I 5270A and 5406A) 
are seen as well as Mg lb 5 172 A which is partly detected at 
the blue edge of the atmospheric absorption feature (the A 
band) . 

For B 2045+265, although signal-to-noise ratio of the 
observed spectrum is worse than that of HST 14113+5211 
due to both fainter apparent brightness and shorter expo- 
sure time achieved for the galaxy, we have clearly detected 
Ca II H and K absorption lines and G band feature (Fig 
We note, however, that these features are less promi- 
nent than those of HST 14113+5211. We also note that 
B 2045+265 shows another possible prominent absorption 
feature at just red edge of the A band feature. We checked 
this possibility by comparing the A band features of both 
galaxies in the observed frame (at 7660 A), and confirmed 
that B 2045+265 shows an excess absorption over expected 
A-band feature. This feature, at 4100A in the rest frame, 
can be identified as H<5 absorption, although another Balmer 
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Figure 2. Same as figure HI but for B 2045+265. The sky spec- 
trum is flux-scaled by a same amount as used for figure Ex- 
pected contribution of the galaxy absorption over the atmospheric 
A band absorption near 4100A is shown in red, with a vertical 
offset of —0.3, overlaying the Gaussian-broadened (a = 213 km 
s — 1 ) template spectrum of HD 136202. See a main text for details 
of this feature. 

absorption line, H7 near the G band, was not detected. 
Since H7 emission is ~ 1.8 times brighter than H5 emis- 
sion under the "case B" photoionization, which is typical 
for star-forming regions, difference in the observed prop- 
erties of Balmer lines can be naturally understood. Note 
that more higher-order Balmer absorptions were identified 
by Fassnacht et al. (1999). Prominent [Oil] emission is also 
detected, suggesting a star-forming activity in this galaxy. 
Fassnacht et al. (1999) classified the spectrum as the Sa 
type, and are best represented by giant stars of late F 
type. The redshift of the lensing galaxy is measured to be 
2 = 0.8682 ± 0.0001 by [Oil] emission, which is consistent 
with that measured by prominent absorption lines, although 
it is slightly larger than the value measured by Fassnacht et 
al. (1999) (z = 0.8673 ± 0.0005). 

2.4 Fourier cross-correlation analysis 

We basically followed the procedure of Ohyama et al. (2002) , 
and the Fourier cross-correlation method (Tonry & Davis 
1979) was used to measure the line-of-sight velocity disper- 
sion of the lensing galaxies with the FXCOR task imple- 
mented in IRAF (see also Falco et al. 1997). Here we outline 
our procedure applied for both galaxies. Firstly, velocity- 
resolution matched stellar spectra were created from the 
"Coude feed spectral library" spectra (Leitherer et al. 1996), 
whose wavelength resolutions are 1.8A FWHM for both red 
and blue spectra in the library, to match the resolution of 
the observed galaxy spectra, and were used as "templates" 
which are to be fitted with the observed galaxy spectra. Sec- 
ondly, calibration curves, which relate the velocity dispersion 
of Gaussian broadening function applied to the templates 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Properties of template stars: (1) star name, (2) spectral 
type, (3) effective temperature (K), (4) surface gravity in log and 
(5) metal abundance in units of solar metallicity. 



StarW 


TypeW T eff (3) 


log gW 


[Fc/H]( 5 ) 


for HST 14113+5211 


HD 206453 


G8III 5092 


2.14 


-0.42 


HD 1918 


G9III 4865 


2.01 


-0.47 


HD 126778 


K0III 4847 


2.34 


-0.53 


for B 2045+265 


BD+11 2998 


F8III 5425 


2.30 


-1.17 


HD 136202 


F8III-IV 6030 


3.89 


-0.07 



to the width of the cross-correlation functions (CCF) peak, 
were created. Finally, curves were used to find the velocity 
dispersion of the galaxy from the CCF peak width calculated 
between the observed galaxy spectra and the templates. FX- 
COR parameters were kept unchanged from that used by 
Ohyama et al. (2002) except for template stars and fitting 
wavelength regions to match the spectral properties of the 
program galaxies in this work. 

2.4.1 HST 14113+5211 

For HST 14113+5211, we choose two separate wavelength 
regions (4000 - 4450 A and 4780 - 5170A, hereafter blue 
and red fitting regions, respectively) to be used in the FX- 
COR, which include most of the prominent absorption fea- 
tures, such as G band, Fe I, and H/3, while avoiding atmo- 
spheric absorption features (A and B bands). Fe I features 
at 5250 - 5410A are not included in the analyses because 
they are detected at redder part of the A band and avail- 
able wavelength region around them is not wide enough. The 
Ca II H and K features were also excluded from the fitting 
regions, since the calculation with these features might give 
problematic result because of their intrinsically wider line 
width (Tonry 1998). We choose late type giant stars with 
sub-solar metallicity as spectral templates, since we found 
during the course of the spectrum fitting that stars with 
nearly solar metallicity show deeper Fe I absorptions than 
the observed ones. We choose three such stars, HD 206453, 
HD 1918, and HD 126778 (see Table □ for their detailed 
properties) 2 . For each star, blue (3820 — 4500A) and red 
(4780 - 5450A) spectra are available in the spectral library, 
and each of them was used separately in the FXCOR for 
fitting blue (4000 - 4450A) and red (4780 - 5170A) fitting 
regions, respectively. Red and blue template spectra were 
processed separately in both the continuum normalization 
and velocity resolution matching procedures, and then each 
template was used separately in the FXCOR to make two 
calibration curves for a set of one template star and the 
galaxy. 

We made FXCOR runs twelve times for each template 
(three stars) , fitting region (two regions per star) , and CCF 

2 Metallicity of the selected template stars ([Fe/H]= —0.4 0.5) 

is by no means the best estimated value for the galaxy. They are 
arbitrarily chosen from a small subsample of low-metal stars in 
the spectral library. 
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Table 2. Results of Fourier cross-correlation analyses for HST 14113+5211. (1) The blue and red fitting regions are 4000 — 4450A, and 
4780 - 5170A respectively. (2) R denotes Tonry-Davis R value. 



Star 


CCF fitting width 


a (km s 1 ] 


bluet 1 ) 


redM 
a (km s — ) 


RM 


HD 206453 


30 


173 


16.6 


162 


13.7 




35 


172 


16.6 


159 


13.6 


HD 1918 


30 


175 


15.5 


179 


14.0 




35 


172 


16.6 


170 


14.1 


HD 126778 


30 


191 


17.8 


170 


14.7 




35 


190 


18.0 


168 


14.8 


R- weighted average a: 
for each fitting region 
for both fitting regions 




179±9 km s" 


-1 

174±10 km s" 1 


168±7 km s" 1 





peak width measurement lag (two lags, 30 and 35, per set 
of template and fitting region; see Ohyama et al. 2002 for 
details of the parameter) to measure velocity dispersion of 
HST 14113+5211. We found that all cross-correlations be- 
tween the galaxy and templates gave good Tonry-Davis R 
values (R = 14—18; Tonry and Davis 1979), suggesting that 
templates adopted for HST 14113+5211 could represent the 
galaxy spectrum well. We found that error of the galaxy 
velocity dispersion, originating from the error in measuring 
the CCF peak widths, was typically ~ 10 km s -1 (la). Av- 
eraging all the FXCOR results with a weight of Rs, we ob- 
tained a — 179 ± 9 km s _1 and 168 ± 7 km s _1 in blue and 
red fitting regions, respectively. Here, errors represent la 
scatters of the velocity dispersion values among all FXCOR 
runs. Another source of error, the velocity resolution match- 
ing error, estimated from the error in measuring the spectral 
resolution of the galaxy spectrum, could amount to ~ 14 km 
s^ 1 and ~ 7 km s _1 of galaxy velocity dispersion error (la) 
for blue and red fitting regions, respectively. Considering all 
sources of uncertainties (errors in velocity-resolution match- 
ing, CCF peak width measurement, and scattering of the 
velocity dispersion values among FXCOR runs), the overall 
la uncertainties are estimated to be 19 (V14 2 + 10 2 + 9 2 ) 
and 13 (V7 2 + 10 2 + 7 2 ) km s _1 for blue and red fitting 
regions, respectively. Note that results measured in two fit- 
ting regions are consistent to each other within expected 
errors. Therefore, averaging all twelve FXCOR results made 
in blue and red fitting regions, we adopt a — 174 ± 20 km 
s" 1 as the best estimated line-of-sight velocity dispersion of 
HST 14113+5211. Here, an error is estimated, based on ve- 
locity matching error of 14 km s _1 (larger value out of two 
fitting regions), CCF fitting error of 10 km s -1 , and scatter- 
ing among all FXCOR results (10 km s , see Table|2J, to 
be 20 (V14 2 + 10 2 + 10 2 ) km s" 1 . 

Fig [T] compares the spectrum of HST 14113+5211 with 
the best template, HD 126778, giving the best R values for 
both blue and red fitting regions, after Gaussian-broadened 
with velocity dispersion of a = 174 km s~ . We also show 
residual spectrum (galaxy — Gaussian-broadened template) 
along with the sky spectrum. All the spectra are shown over 
wider wavelength region, covering other spectral features not 
used for the FXCOR analyses, to show the overall quality 
of the spectral fitting. One may find that the broadened 
template gives rather good representation of the observed 
galaxy spectrum over entire wavelength region, except for 



a large-scale residual pattern in the blue continuum. The 
residual feature likely comes from errors in fitting the con- 
tinuum spectra of both the galaxy and the templates, used 
for continuum normalization. Since the continuum fittings 
were made over the wavelength regions where spectra show 
rather large-amplitude change in shape, an artificial large- 
scale spectral variation could be easily introduced into the 
continuum-normalized spectra. Note, however, that this er- 
ror does not affect our cross-correlation results, since such 
a low-frequency variation in spectral shape is filtered out 
in the FXCOR before calculating the CCF (see Ohyama et 
al. 2002 for more details on the FXCOR parameters). In 
concluding, after considering the fitting quality as discussed 
above, we finally adopt a = 174 ± 20 km s _1 as the best 
estimated velocity dispersion of HST 14113+5211 in the fol- 
lowing lens analysis. 



2.4.2 B 2045+ 265 

For B 2045+265, we followed the same procedure used 
for HST 14113+5211. However, because of relatively poor 
quality of B 2045+265 spectrum, only a narrow wave- 
length region, comprising two smaller wavelength regions of 
4280-4325A and 4355-4430A, was chosen for the FXCOR 
analyses. The spectral regions were carefully chosen to in- 
clude G band while avoiding both the strong sky emission 
regions and H7, and these two smaller regions were fitted at 
a same time in FXCOR. Here, the H7 was excluded in the 
fitting since H7 emission may overlay underlying H7 absorp- 
tion. We choose two late F giant stars, BD+11 2998 and HD 
136202, as spectral templates for the galaxy from the same 
spectral library (see Tabled for their detailed properties) 3 . 
We ran FXCOR four times (for one fitting region, two tem- 
plate stars, and with two CCF peak width measurement 
lags), and found that R values are rather large, 30 — 40, 
for all FXCOR runs. The larger R values probably result 
from the fact that we have only one prominent spectral fea- 
ture (G band) in a narrow wavelength region for fitting. The 
best estimated velocity dispersion is a = 213 km s , with 
an la error of ±11 km s _1 originating from scattering of 

3 Again, as in the case of the analyses of HST 14113+5211, mctal- 
licity of the selected template stars is by no means the best esti- 
mated value for the galaxy. 
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Table 3. Same as TableEJbut for B 2045+265. Fitting region is 4280 - 4325A plus 4355 - 4430A. 



5 



Star 


CCF fitting width 


a (km s 1 ) 


R 


BD+11 2998 


30 


219 


30.9 




35 


234 


30.5 


HD 136202 


30 


210 


39.9 




35 


213 


39.7 


i?-weighted average <r: 




213+11 km s" 1 





the a values among all FXCOR runs (see Table 0J. The ve- 
locity resolution-matching error could amount to ~ 5 km 
s -1 of galaxy velocity dispersion error (la), which is smaller 
than in the case of HST 14113+5211 due to higher redshift 
of B 2045+265. We found that error of the galaxy velocity 
dispersion, originating from the error in measuring the CCF 
peak widths, was typically ~ 19 km s -1 (la), which is larger 
than the value for HST 14113+5211 because of relatively 
poor CCF determination resulting from lower quality spec- 
trum and narrower wavelength fitting regions available for 
B 2045+265. Considering all sources of uncertainties (errors 
in velocity-resolution matching, CCF peak width measure- 
ment, and scattering of the velocity dispersion values among 
FXCOR runs) , the overall la uncertainty is estimated to be 
23 (V5 2 + 192 + ll 2 ) km s -1 . 

Fig [5] compares the observed galaxy spectrum with the 
Gaussian-broadened template spectrum of the best tem- 
plate, HD 136202, after Gaussian-broadened with the ve- 
locity dispersion of a — 213 km s -1 , as in the same way for 
HST 14113+5211. One may find that the broadened tem- 
plate gives rather good representation of the galaxy spec- 
trum over entire wavelength region, including US, except for 
H7 and Ca II H and K features. We expect that the residual 
emission of H7 represents the real emission, overlaying on 
the absorption line, as expected from another emission line 
of [Oil]. On the other hand, it seems difficult to attribute 
the shallower-than-observed Ca II absorptions in the broad- 
ened template spectrum to overlaying emissions. Although 
putative He emission, whose wavelength is almost coincident 
with that of Ca II H absorption, might explain a part of the 
residual feature, there seems no emission lines which could 
be coincident with Ca II K feature neither in the spectra 
of the lens galaxy nor the lensed QSO. Therefore, the most 
simple and plausible explanation for the shallower observed 
Ca II absorption is that the lensing galaxy contains not only 
stars of late F type but also of even earlier types. Although 
neglecting the contribution of the earlier population in the 
FXCOR analyses should affect the velocity dispersion mea- 
surement, we expect that our result, made assuming a pure 
stellar population of late F type as spectral templates, is 
still reasonably reliable and is useful for our analyses on the 
mass model and the stellar dynamics in the following sec- 
tions. This is because (1) stars of earlier types show much 
shallow G band feature, and they can not fit the observed 
spectrum alone, indicating that late F stars contribute most 
to the observed spectrum in flux, and (2) our measurement 
almost relies on the spectral fit of the G band feature and, 
hence, the velocity dispersion measured in our analyses lies 
larger weight on the later type stars, which contribute much 
more on the mass of the lensing galaxy than the case of 
stars of earlier types. Therefore, we adopt a — 213 ± 23 km 



s 1 as the best estimated line-of-sight velocity dispersion of 
B 2045+265 in the following lens analysis. 



3 FUNDAMENTAL PLANE 

Before going into a joint gravitational lensing and stellar 
dynamics, let us examine the Fundamental Plane (FP) for 
early type galaxies, which provides us with an alternative in- 
sight into the value of the mass-to-light ratio M*/Lb- The 
FP is defined by the effective radius (Re kpc), central ve- 
locity dispersion (a c km s -1 ), and mean surface brightness 
inside R e (SB e mag arcsec -2 ): 

log R e = a FP log a c + j3 FP SB e + 7fp , (1) 

where off = 1.25, Pfp = 0.32, and jfp = —8.895 — log hso 
(hso = Ho/50 km s -1 Mpc -1 ) in the B band (Bender et 
al. 1998). The central velocity dispersion a c is taken as the 
velocity dispersion measured inside R B /8 (TK04). 

For the lensing galaxy of HST 14113+5211, Fischer, 
Schade, & Barrientos (1998, hereafter F98) obtained the 
rest-frame central surface brightness of the lensing galaxy 
as fio,B(AB) — 13.57 mag arcsec -2 . Then, as the surface 
brightness at R e = 0."61 ± 0."03 (F98) is fi e , B (AB) = 
fi , B (AB) + 8.33 and B(Vega) = B(AB) + 0.12, we obtain 
the surface brightness at R e as \ie,B = 22.02 mag arcsec -2 . 
The mean surface brightness within R e is then estimated 
as SB e — 20.56 mag arcsec -2 after correction for Galactic 
extinction E(B — V) = 0.016 mag (Schlegel, Finkbeiner, & 
Davis 1998) and an Ry = 3.1 extinction curve. Alterna- 
tively, van de Ven, van Dokkum, & Franx (2003, hereafter 
vdV03) derived /i e ,B = 21.70 mag arcsec -2 at the intermedi- 
ate effective radius R e of 0."48 (Kochanek et al. 2000), giv- 
ing SB e = fi e ,B — 1.393 ~ 20.56 mag arcsec -2 . As the small 
difference between these values results in only the difference 
of 0.2 to 0.3 in the mass-to-light ratio obtained below, we 
adopt the results of F98 in what follows. 

For the lensing galaxy of B 2045+265, we adopt the 

values given in vdV03, i.e., R e = 0."38 + °„^ and SB e = 

/i e ,s - 1.393 ~ 19.11 ± 0.44 mag arcsec -2 . 

To derive a fiducial correction factor from the measured 
a in a specific aperture to a c in equation Q, we utilize a 
fiducial model for internal mass distribution and stellar dy- 
namics as explained in Briefly saying, we adopt a single 
mass component represented by a power-law density pro- 
file with an index of 7' [see equation (5)] and the veloc- 
ity dispersion of stars represented by a constant anisotropy 
parameter (3 [see equation (6)]. For HST 14113+5211, we 
obtain a c — 1.15<r based on the most likely parameters of 
(7',/?) = (1.93,0.14) (see © thereby yielding a c = 200.7 
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Figure 3. Lens configuration of the HST 14113+5211 compared 
with SIE lens with an external shear model. Multiple image po- 
sitions and source position of the best-fitting model are shown 
by crosses and plus, respectively. Open circles show the observed 
image positions, where the radii of the circles correspond to 1-cr 
observational uncertainties. The lensing galaxy position is shown 
by filled circle and is taken as the coordinate origin. 

km s _1 . For B 2045+265, we obtain a c — 0.99a based on 
the most likely parameters of (y',/3) = (1.66,0.10), thereby 
yielding cr c = 210.3 km s . 

Inserting into equation ifT), we obtain 7pp = —8.87 
for the lensing galaxy of HST 14113+5211 at z L = 0.464 
and 7fp = —8.52 for the lensing galaxy of B 2045+265 at 
zl = 0.87. Comparison with 7pp at z — allows us to derive 
the difference in 7fp, A7PP, and then Alog(M*/Lp) = 
—OAAjfp/Pfp provided «pp and (3fp are constant. We 
then obtain Alog(M,/Lp) = —0.172 for the lensing galaxy 
of HST 14113+5211 and -0.612 for the lensing galaxy of 
B 2045+265. The mass-to-light ratio at z > is given as, 

log(A/,/Lp) z = log(M,/L s ) + A \og(M/L B ) (2) 

where (M*/L B )o = 7.3 ± 2.1/i 65 M Q /L B ,e (Gerhard et al. 
2001). We obtain (M*/L B )* = 4.9 ± 1.4 M Q /L B ,e for the 
lensing galaxy of HST 14113+5211 and 1.8 ± 0.5 M /L s , 
for the lensing galaxy of B 2045+265. As will be shown be- 
low, these values of M*/Lp for both systems are in good 
agreement with those obtained from a joint lensing and dy- 
namical analysis (see [0. 

4 MASS MODEL AND STELLAR DYNAMICS 
4.1 Lens models 

Following TK04, we model the lensing galaxy as a singu- 
lar isothermal ellipsoid (SIE: Kormann, Schneider & Bartel- 
mann 1994). Note that the Einstein radius (Re) and mass 
enclosed by the Einstein radius (Me), both of which are 
the quantities required in the dynamical model of the lens- 
ing galaxy in the following sections, are very insensitive to 
the assumed mass model (Kochanek 1991; Koopmans & 
Treu 2004). We also allow for a constant external shear. 
The observed position of a lensing galaxy is taken by a 
lens position and we do not treat the lens position as a 
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Figure 4. Same as Fig. EJ but for B 2045+265. 



Table 4. Best-fitting model parameters and characteristic values 
of the singular isothermal ellipsoid with an external shear lens 



Parameter 


HST 14113+5211 


B 2045+265 


b[ (arcsec) 


0.84 


1.11 


H 


0.68 


0.68 


61 (deg) 


38 


-69 


7oxt 


0.26 


0.05 


cx t (deg) 


-35 


18 


j3 x (arcsec) 


-0.005 


-0.39 


f} y (arcsec) 


0.11 


0.12 


cr SIE (km/s) 


202 


397 


Re (kpc) 


5.29 


9.19 


M(< R E )(M Q ) 


1.58 x 10 11 


1.06 x 10 12 



free parameter. Therefore, our lens model has five parame- 
ters (we follow the notations in TK04); the lensing strength 
(6; = 47t((jsie/c) 2 U; s /D s , where asm is the one-dimensional 
velocity dispersion of SIE lens), axis ratio (qi), position an- 
gle of the lens (Qi), strength of the external shear (7 ex t), and 
its orientation (# ex t). Also we treat the source position (f3 x , 
/3 y ) as a free parameter. 

We search for a set of model parameters that best re- 
produces the observed lens configuration. The best-fitting 
model parameters are summarized in Table 21 Fig s CD 
and 0] compare the observed image positions with the 
model images. The usual \ 2 values are 0.71 and 14 for 
HST14113+5211 and B 2045+265, respectively (for l-a ob- 
servational uncertainties of 0.03 arcsec) and the number of 
degrees of freedom is Ndof = 1 (8 constraints and 7 param- 
eters). As is evidently shown in FigEH for HST14113+5211, 
the model reproduces the lens configuration very well. The 
Einstein radius and mass enclosed by the Einstein radius are 
found to be i?p=5.29kpc and M(< R E ) = 1.58 x 10 n M Q 
(for the adopted cosmological parameters), respectively. On 
the other hand, for B 2045+265, although the model nicely 
reproduces the cusp-lensing configuration (Schneider Ehlers 
& Falco 1992), the best- fitting model positions slightly de- 
viate from the observed positions. Since the enclosed mass 
is not very sensitive to a detail lens model but is primarily 
determined by the image separation, the best-fitting model 
may give a good estimate of the the enclosed mass. The Ein- 
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stein radius and mass enclosed are found to be 7?g=9.19kpc 
and M(< Re) = 1.06 x 1O 12 M , respectively, which we 
adopt in the following joint lensing and stellar dynamics 
analysis. 



4.2 Mass model and stellar dynamics 

Let the luminous- mass density and dark-matter density pro- 
files be pium(r) and pdm(t), respectively, provided that 
these are spherically symmetric. Two alternative mass mod- 
els will be employed, following TK04. 

First, we consider a two-component mass model, where 
the luminous mass and dark matter distribute differently: 



pium(r) — 



Pdm(t) = 



2nr(r + r„) 3 
Pdm,otI 



r 2\(3- 7 )/2 



(3) 



(4) 



where M, is the total stellar mass and r* denotes the scale 
length for the luminous matter. The profile piumir) corre- 
sponds to a Hernquist model (Hernquist 1990), reproducing 
the _R 1//4 surface brightness profile with an effective radius 
of R e = 1.8153r*. For M», we use the B-band total lumi- 
nosity of the lensing galaxy, Lb , based on the mass-to-light 
ratio of the luminous matter, M*/Lb, as a model param- 
eter. Our estimation of Lb is given in Appendix, yielding 
Lb/L b ,q = 2.6 x 10 10 and 1.3 x 10 11 (h = 0.65) for the 
lensing galaxies of HST 14113+521 and B 2045+265, re- 
spectively. The dark-matter density profile, pdm(v), is de- 
termined by the scale length (r b ), density scale (pdm,o), and 
inner slope (7). The combination of all of these parameters 
will be constrained by the results of the lens fitting, Re and 
M(< Re), as obtained in §4.1. 

Second, we consider a single component with power-law 
mass model: 



Ptot(r) oc r 



(5) 



where 7' denotes an effective slope. 

The velocity distribution of the stars is also assumed 
to be spherically symmetric, such that the velocity disper- 
sions in spherical coordinates (cr r ,crg,a^) satisfy ag — a^, 
(e.g., Binney & Tremaine 1987). We then use the parameter 
f3(r) = 1—o-g /crj? describing the degree of velocity anisotropy. 
We employ the following two models for f3(r), the Osipkov- 
Merritt model with a parameter rt (referred to as Model A) 
and constant anisotropy model with a parameter /3 (Model 
B). 



0(r) 



const. (= /3), 



for Model A 
for Model B 



(6) 



Relevant parameters are n and (3. 



5 RESULTS OF A JOINT LENSING AND 
DYNAMICAL ANALYSIS 

Based on the models given in the previous section, we 
search for the best mass model for reproducing the ob- 
served velocity dispersion a using the joint lensing and 
dynamical analysis and also combining with the FP con- 
straints. We especially focus on the best values of M*/Lb 



M.O 



M.O 




Figure 5. Likelihood contours for a joint lensing and dynam- 
ical analysis of HST 14113+5211 based on the two-component 
mass model. Contours show 68 % (la), 95 %, 99 %, and 99.9 % 
confidence levels. Dashed lines are derived from the current joint 
analysis, while solid lines combine the additional constraints from 
the FP plane. Lower panel: Model A with t\ = R e and = Re- 
Upper panel: Model B with /3 = and r\y = Re- 



for a two-component model, which are to be compared 
with those from the FP constraints. Also, the inner slope 
of the dark-matter halos, 7, are to be compared with the 
prediction of cosmological N-body simulations [e.g., 7 = 
1 - 1.5 (Navarro, Frenk, & White 1995; Moore et al. 1998; 
Fukushige & Makino 2003 and references therein)]. For a 
single-component model, the derived slope 7' can be used 
to assess the isothermality of the lens density profile as usu- 
ally adopted in other lensing work. 
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Figure 6. Likelihood contours for a joint lensing and dynamical 
analysis of HST 14113+5211 based on the single component with 
power-law mass model. Contours show 68 % (lc), 95 %, 99 %, and 
99.9 % confidence levels. Lower and upper panels are for Model 
A and Model B, respectively. 



5.1 HST 14113+5211 

First, we examine a two-component mass model. Figure 
shows 68 % (la), 95 %, 99 %, and 99.9 % confidence levels of 
the likelihood contours in the plane of 7 and M*/Lb when 
we adopt the velocity dispersion of a = 174 ± 20 km s _1 . 
Dashed lines are derived from the joint lensing and dynam- 
ical analysis, while solid lines combine the additional con- 
straints from the FP plane. For the lower panel (Model A) 
we set n = R e and rt = Re, and for the upper panel (Model 
B) we set /3 = and r& = Re- The likely 7 or M»/Lb us- 
ing these parameter sets will be utilized as a characteristic 



case in what follows, as the likelihood contours are found 
to remain basically the same in other parameter settings 
(TK04). 

From the results based on the current joint analysis 
(dashed lines), the likely mass-to- light ratio appears to be 
in the range of 4 < M*/Lb < 6 Mq/ Lb,q'- after marginal- 
izing over 7, we obtain the most likely values of M*/Lb 
(M q /L b ,q), yielding 4.4t2: 1 for Model A and 5.6±l% for 
Model B. As is evident, these values of M*/Lb are virtually 
consistent with those obtained from the FP in §3 (4.9 ±1.4 
Mq/Lb,q)- We also calculate the most likely fraction of the 
dark matter projected inside Re (— lAR e ), which is de- 
noted as /dm- Adopting the isotropic velocity model as a 
representative case, we obtain /dm = 0. 471^2!, suggesting 
that about a half of the total mass derived from the lens 
model constitutes dark matter. 

Turn next to the result from the joint analysis combined 
with the FP constraints (plotted by the solid lines). Since 
the above joint analysis gives M*/Lb being in a good agree- 
ment with one from the FP, the additional constraint from 
the FP makes only a minor change in the preferred M*/Ls 
value. It is found that marginalized constraints on M*/Lb 
(Mq/Lb,q) are 4.8+^ for Model A and 5.1±^ for Model B. 
The most likely fraction of the dark matter projected inside 
Re (for the isotropic velocity case) reads f DM = 0.52±° 
The additional constraint from the FP improves the limit 
on 7. We obtained, after marginalizing over M,/Lb, a con- 
straint on 7 to be 7 < 1.6 (7 < 1.8) at 1-cr for Model A 
(Model B). 

Second, we examine a single-component mass model to 
set limits on the slope 7'. Figure|S|shows 68 % (lcr), 95 %, 99 
%, and 99.9 % confidence levels of the likelihood contours. 
We set Tb = 3Re in this diagram, but the limits on 7' are 
found to be little affected as long as rt > 3Re- We esti- 
mate the most likely values of 7' for n = R e (Model A) and 
(3 = (Model B), yielding 1.87±° ° l and 1.94±g;^, respec- 
tively. Thus, the total density profile is well approximated 
as p t ot(r) oc r~ 2 . 



5.2 B 2045+265 

For the lensing galaxy of B 2045+265, we show the results 
in Figure □ and Q3 which are to be compared with Figure 
and© respectively, obtained for HST 14113+5211. 

For a two-component mass model (Figure 0, the likely 
mass-to-light ratio appears to be in the range of 0.7 < 

we 
for 

Model A (Model B) as M»/L B = 0.8toi (l.Z±$5) without 



M*/Lb J; 2 Mq /Lb,®'- after marginalizing over 7, 
obtain the most likely values of M*/Lb {Mq/ Lb,q) 



the FP constraints. As is evident, these values of M*/Lb 
are also consistent with those obtained from the FP in §3 
(1.8 ± 0.5 Mq/Lb,q)- Combining the FP constraint and af- 
ter marginalizing over M*/Lb we obtain 7 < 0.5 (7 < 0.8) 
at 1 a for Model A (Model B). Thus the models prefer flatter 
dark matter inner slope than the case of HST 14113+5211. 

We also obtain the most likely fraction of the dark mat- 
ter projected inside Re (— 2.9-R e ) for the isotropic velocity 
case, yielding }dm = 0.891q'q4 (without the FP constraints) 
and 0.86±°;°s (with the FP constraints). Thus, in contrast to 
the case of HST 14113+5211, the total mass inside Re de- 
rived from the lens model is totally dominated by dark mat- 
ter; some of this unseen mass component may be provided 
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Figure 7. The same as FigEJbut for B 2045+265. 
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Figure 8. The same as Fig|H]but for B 2045+265. 



by the group of galaxies which the lensing galaxy belongs 
to. 

For a single-component mass model (Figure [HJ, we ob- 
tain the most likely values of 7' for n — R e (Model A) and 
/3 = (Model B), yielding 1.58±g;g| and 1.66±g;g|, respec- 
tively. Thus, the model prediction for the slope of the total 
density profile is systematically flatter than isothermal. 



6 DISCUSSION AND CONCLUDING 
REMARKS 

As explored in this work and also in TK04, the mass-to- 
light ratio for both of our targets based on the joint lensing 
and dynamical analysis is virtually consistent with that ob- 



tained from the FP constraint; the mass-to-light ratios at 
zl = 0.464 and 0.868 are systematically smaller than the 
current average value of (M»/Lb)o = 7.3±2.1/i65 Mq/Lb,q 
(Gerhard et al. 2001), implying the aging of stellar popula- 
tions from these redshifts to the present day. For comparison 
with other sample lenses for which the similar analysis has 
been employed, we plot, in Figure^ the redshift evolution of 
the stellar mass-to-light ratio for both of our targets (open 
circles) and the TK04 sample (solid circles). It is found that 
the current sample lenses follow the general redshift evo- 
lution of Mt/Ls, as guided by the solid and dashed lines 
showing the average evolution of M*/Lb derived from the 
FP of early-type lensing galaxies by TK04 and Rusin et al. 
(2003), respectively. 

The inner slope of the dark-matter halos, 7, ranges from 
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Figure 9. Redshift evolution of the stellar mass-to-light ratio. 
The open circles correspond to HST 14113+5211 and B 2045+265 
derived from the current lensing model (for an isotropic velocity) 
without the FP constraints, whereas the solid circles are taken 
from TK04 with the same models. The solid and dashed lines show 
the average evolution of M* /Lb derived from the FP of early- type 
lensing galaxies by TK04 and Rusin et al. (2003), respectively. 



0.0 to 1.8. More specifically, the most likely values of 7 (with 
the FP constraints) for HST 14113+5211 are < 7 < 1.6 
(0 < 7 < 1.8) at 1 a for Model A (Model B), whereas for 
B 2045+265, we obtain < 7 < 0.5 (0 < 7 < 0.8) at 1 a 
for Model A (Model B). These values are consistent with 
the density slopes of dark halos predicted by cosmological 
N-body simulations, showing 7 of 1.0 to 1.5 (e.g., Navarro, 
Frenk, and White 1995; Moore et al. 1998; Fukushige & 
Makino 2003 and references therein), although the present- 
day inner slope of the dark-matter halos is affected by the 
effect of baryonic infall to some extent; the effect seems to 
be not so significant as obtained from other lens samples 
(TK04). It is also drawn that the inner slope of the dark- 
matter halos is systematically flatter than an isothermal pro- 
file. 

The index for a single power-law model (7') is almost 
2 for HST 14113+5211, thereby suggesting that a singular 
isothermal model is a good fit to this lensing galaxy. For 
B 2045+265, 7' ranges from 1.5 to 1.7, which is system- 
atically flatter than isothermal (Figure [HJ. This flat slope is 
related to the very small velocity dispersion of stars a — 213 
km s _1 compared with asm = 397 km s - ; a flatter slope of 
a dark halo than SIE is needed to reduce the radial gravita- 
tional force and thus the velocity dispersion of stars. Alter- 
natively, some part of M(< Re) for B 2045+265 is provided 
by the group of galaxies, thereby causing a large value of 
asm. This effect may partly contribute to a non-negligible 
scatter of 7' from an isothermal index (2) as is also reported 
in other lens samples (TK04). 
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APPENDIX A: ESTIMATION OF THE LENS 
LUMINOSITY IN THE B BAND 

For the lensing galaxy of HST 14113+5211, F98 reported the 
F702W(AB) magnitude as 20.78±0.05 mag. Based on the K- 
correction using Coleman, Wu, & Weedman (1980) and the 
assumption of no galaxy evolution, F98 derived the absolute 
B magnitude as M B (AB) = -19.32 + 5 log ft for (n ,A ) = 
(1,0) and zl = 0.46. Then, we transform this magnitude 
into Vega-based B magnitude [M s (Vega) = M B (AB)+0.12, 
Schmidt, Schneider, & Gunn 1995] and consider the revised 
lens redshift of zl = 0.464 (Lubin et al. 2000) and cosmolog- 
ical parameters of (Oo, Ao) = (0.3, 0.7) (Spergel et al. 2003). 
We thus obtain Mb = —19.60+5 log h, giving the luminosity 
of Lb/Lb,q = 2.6 x 10 10 (h = 0.65) for the lensing galaxy 
of HST 14113+5211. 

For the lensing galaxy of B 2045+265, Fassnacht et al. 
(1999) reported the various infrared magnitudes in a l."9 
diameter aperture (corresponding to the size of the Ein- 
stein ring in their lens model). Based on the K-correction 
for an Sa galaxy as the lens appears to show its typical 
spectrum, they arrived at the rest-frame B-band luminosity 
of 2.36 x 10 10 /i" 2 L o in this aperture for (f2 ,A ) = (1,0). 
To derive the total B-band luminosity, we adopt the work of 
Rusin et al. (2003), where they obtained the intermediate- 
axis effective radius (R e ) determined by fitting the observed 
brightness distribution to a de Vaucouleurs profile. Using 
their R e of 0."38 and (fi , A ) = (0.3, 0.7), we obtain the lu- 
minosity of Lb/Lb,q = 1-3 x 10 11 (h = 0.65) for the lensing 
galaxy of B 2045+265. 
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